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Abstract 

The concept of p-adic quincunx Haar MRA was introduced and stud- 
ied in [6]. In contrast to the real setting, infinitely many different wavelet 
bases are generated by a p-adic MRA. We give an explicit description for all 
wavelet functions corresponding to the quincunx Haar MRA. Each one gen- 
erates an orthogonal basis, one of them was presented in [6]. A connection 
between quincunx Haar bases and two-dimensional separable Haar MRA is 
also found. 
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1 Introduction 

In 2002, S. V. Kozyrev [13] found ap-adic wavelet basis for L 2 (Q p ) which is 
an analog of the real Haar basis. V. M. Shelkovich and one of the authors de- 
veloped an MRA approach to p-adic wavelets [TTj. The scheme was realized 
to construct the p-adic Haar MRA. In contrast to the real setting, it appears 
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that there exist infinitly many different orthonormal wavelet bases with the 
minimal number of generating wavelet functions in the same Haar MRA. 
All these bases were described explicitly in [17] for p = 2 and in [11] for 
arbitrary p, one of these bases (for each p) coincides with Kozyrev's wavelet 
basis. A wide class of orthogonal scaling functions generating an MRA was 
described in [10] by A. Khrennikov, V. Shelkovuch and one of the authors. 
However it has been shown in [1] by S. Evdokimov and both authors that 
all these scaling functions lead to the same Haar MRA and that there exist 
no other orthogonal test scaling functions generating an MRA, except for 
those described in |10j . Thus all univariate orthogonal MRA-based wavelet 
bases with a minimal possible number of wavelet functions are described. It 
is not known if other p-adic orthogonal wavelet bases exist. Although, for 
the construction p-adic wavelets, there is another approach based on p-adic 
wavelet set theory introduced by J. J. Benedetto and R. L. Benedetto [3], 
[2] , we did not see up to now orthogonal p-adic wavelet bases which are not 
generated by the Haar MRA. 

The simplest way to construct a multivariate wavelet basis is the follow- 
ing. Given an univariate MRA with a scaling function 93, using the method 
suggested by Y. Meyer p3] (see, e.g. [151 Sec 2.1]), one can easily construct a 
c?-dimentional separable MRA, where the scaling function is <£>(g> • • • <S>ip. Hip 
is a wavelet function in the univariate MRA, then the functions fx (g> ■ ■ ■ (g) fa, 
where each fk equals either ip or ip and not all k = 1, . . . , d, are equal 
to <y9, form a set of wavelet functions generating a multivariate wavelet basis 
(i.e. the basis consists of the dilations and shifts of the wavelet functions). 
This approach was realized for p-adics in |17] . Next we are going to discuss 
only the case p = d = 2. The separable wavelet basis is generated by three 
wavalet functions in this case. E. King and one of the authors studied two- 
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dimensional Haar MRA with quincunx dilation in [6j . The basis is generated 
by only one wavelet function in this case. In the present paper we describe 
all wavelet functions corresponding to the quincunx MRA and generating an 
orthogonal basis. We also discuss how these bases related to the separable 
MRA. 



2 Preliminaries and notations 
2.1 p-adic numbers 

The theory of p-adic numbers is presented in detail in many books (see, 
e.g., [18], [8]). We restrict ourselves to a brief description. 

Let N, Z, M, C be the sets of positive integers, integers, real numbers, 
complex numbers, respectively. The field Q p of p-adic numbers is defined 
as the completion of the field of rational numbers Q with respect to the 
non-Archimedean p-adic norm | • | p . This p-adic norm is defined as follows: 
|0| p = 0; if x ^ 0, x = p 11 ^, where 7 = 7(2;) € Z and the integers m, n 
are not divisible by p, then \x\ p = p~ 7 . The norm | • \ p satisfies the strong 
triangle inequality \x + y\„ < max(|x| p , \y\ p ). The canonical form of any 
p-adic number x 7^ is 

00 

x = p 1 x jp' (1) 
3=0 

where 7 = ^y(x) € Z, Xj £ D p := {0, 1, ... ,p — 1}, xq ^ 0. The fractional 
part {x} p of the number x equals by definition "Y^Iq 1 x jP ' • Thus, {x} p = 
if and only if 7 > 0. We also set {0} p = 0. 

For every prime p, the additive group of Q p is a locally compact abelian 
group which contains the compact open subgroup Z p . The group Z p is 
the set {a G Q p \a\ p < 1}, the unit ball in Q p . Equivalently, 7L p is the 
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subgroup generated by 1. It is well-known (|16j) that since Q p is a locally 
compact abelian group with a compact open subgroup, it has a Haar measure 
normalized so that the measure of Z p is 1. For simplicity, we shall denote 
this measure by dx. 

We also define the set I p = {x G Q p : {x} = x}. I p is a set of coset 
representatives for Q p /Z p . Since 7L p is open, I p is discrete. 

For any d > 1, Q p is a vector space over Q p , i.e., Q p consists of the 
vectors x = (x\, . . . , xd), where Xj G Qp, j = 1, . . . , d. The p-adic norm on 
Qp is 

The Haar measure cix on the field Q p is extended to a Haar measure 
d d x = dx\ • • • dxd on Q p in the standard way. 

Denote by = {x G 
center at the point 0, N G 



1 p 

Denote by B N = {x G Qp : |x|p < p w } the ball of radius p N with the 



Finally, we define the set I p C Qp by I p := I p x • • • x I p . 
2.2 Functions of p-adic variables 

There are two main function theories over the p-adics. One deals with 
functions Q p — > C and the other with Q p — > Q p . We shall only deal with 
the former theory. 

Since a Haar measure exists in Q p , the spaces L 9 (Q p ) can be naturally 
introduced. 

The functions Xp( x ) = e 2m ^ p are additive characters for the field Q p . 
For any set S C Q^, t s denotes the characteristic function of S, i.e. 



l s (x) 



( 1 ; x G S 
; else 



Note that l Zp <g> l Zp = 1 Z 2 . 



2.3 p-Adic MRA 

Definition 1 (|17j). A collection of closed spaces Vj C L 2 (Q P ), j G Z, is 
called a multiresolution analysis (MRA) in L 2 (Q p ) if the following axioms 
hold 

(a) Vj C V j+ i for all j G Z; 

(b) 1J V,- is dense in L 2 (Q p ); 
jez 

(c) fl ^ = {0}; 

j£Z 

(d) /(•) G ^ /(p" 1 -) G for all j G Z; 

(e) there exists a function 99 G Vq, called the scaling function, such that 
the system {(p(- — a), a G I p } is an orthonormal basis for Vq. 

According to the standard scheme (see, e.g., [151 § 1-3]) for the construc- 
tion of MRA-based wavelets, for each j, we define the space Wj (wavelet 
space) as the orthogonal complement of Vj in Vj+i, i.e., 

V j+1 = Vj®Wj, jeZ, (2) 

where Wj -1 Vj, j G Z. It is not difficult to see that 

/ G Wj /(p — 1 -) G Wj+i, for all j G Z (3) 

and _L j 7^ fc. Taking into account axioms (b) and (c), we obtain 

(J)Wj = L 2 (Q P ) (orthogonal direct sum). (4) 

jez 

If now we find functions ip u G Wq, v = 1, ... ,r, such that the system 
{4>v(x — a), a G 7 p , 1/ = 1, . . . , r} is an orthonormal basis for Wq, then, due 
to © and dH, the system {p 7/ ' 2 Vv(p _jf • —a), a G Jp, j G Z, ^ = 1, . . . ,r}, 
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is an orthonormal basis for Z/'((Q)p)« Such, functions i\) v are called wavelet 
functions and the basis is a wavelet basis. 

3 Separable and quincunx ]>adic wavelets bases 
3.1 Separable MR A 

Let {vj }jez > ^ = 1, be one-dimensional MRAs. We introduce 

subspaces Vj, j £ Z, of £ 2 (Qp) by 

d 

^ = (g) vf° = span{F = /i ® • • • ® /„, /„ € Vj W }. (5) 

Let </9 M be a scaling function of z/-th MRA {V^jj. Set 

$ = yjC 1 ) (g) . . . (gj 

It is not difficult to check that the following theorem holds (see |17j). 

Theorem 2. Let {vj u) } j& , v = 1, . . . ,n, be MRAs in L 2 (Q p ). Then the 
subspaces Vj of L 2 (Q^) defined by ([5]) satisfy the following properties: 

(a) Vj C Vj+x for all j 6 Z; 

(b) Dj,zzVj is dense in L 2 (Qp); 

(c) n jeZ Vj = {0}; 

(d) /(•) G Vj <=► /(p- 1 -) G Vj +1 /or aZZ j G Z; 

(e) f/iere exists a function <fr € Vb suc/i £/tat f/ie system {&(x — a), a € i^}, 
forms an orthonormal basis for Vb • 

The collection of spaces Vj, j € Z, discussed in Theorem [2] is called a 
separable multiresolution analysis in L 2 (Qp), the function $ from axiom (e) 
is said to be its scaling function. 
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Following the standard scheme (see, e.g., |X5|, § 2.1]), we define the 
wavelet space Wj as the orthogonal complement of Vj in Vj + i, i.e., 

Wj l). i \' r j€Z. 

Since 

v=l u=l 

eC{l,...,n},e^0 ^Ge 

So, the space Wj is a direct sum of 2 d — 1 subspaces W^e, e C {1, . . . ,n}, 
e/0. 

Let ip^ be a wavelet function, i.e. a function whose shifts (with respect 
to a G 7p) form an orthonormal basis for Wq . It is clear that the shifts 
(with respect to a e /p) of the function 

* e =((g)^ (l/) )((gM M) ), ec{l,...,n}, (6) 

form an orthonormal basis for Wo,e- So, we have 

L 2 (Qj) = ©W, = ©( © W i>e ), 

j&L j&L eC{l,...,n},e^0 

and the functions p~ d i l 2 ^ e {jj ■ +a), e C {1, . . . , n}, e ^ 0, j G Z, a G Ip, 
form an orthonormal basis for L 2 (Qp). A wavelet basis constructed in this 
way is called separable. 

In particular, if p = d = 2, p^ = ip^ = tp, (which implies {Vj }jeZ = 

(2) 

{V^ }jez) and ip(x) = X2{x/2)tz 2 (x) is a wavelet function corresponding 
to the univariate MRA, then there are three wavelet functions X I /1 = (p <S> 
ip, ty 2 = ifi ® 99, ^ 3 = ?/> ® V'- corresponding to the separable MRA. 
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3.2 Quincunx MRA 



Set 

-l 

1 l 



A 



2 2 



1 1 

"2 2 



1 -1 
1 1 

A is the inverse of the well-known quincunx matrix. In |7J (see also \15\ 
§ 2.8]), a Haar multiresolution analysis was presented for L 2 (M 2 ) which used 
dilations by the quincunx. Since | detyl _1 | = 2, there was only one wavelet 
generator. However, the support of the scaling function was a fractal, the 
twin dragon fractal. In contrast, it was proved in |B] that the Haar MRA 
for L 2 (Q 2 ) associated to this matrix corresponds to a scaling function with 
a simple support, namely, Z|. To state this in more detail set 



^ = span{l zi (^--a):ae/ 2 2 }. (7) 

Theorem 3 (©). The sub spaces V9 of L 2 (Q 2 ) ; j G Z d , defined by ^ 
satisfy the following properties: 

(a) V? C V9 +1 for all j G Z, 



(b) UjezV® ^ dense m L 2 (Q 2 2 ), 

(d) ftV®^ f(A-) G V9 +1 for all j G Z, 

(e) there exists a function <f> G V® , such that the system {(/)(■ — a) : a G /|} 
is an orthonormal basis for Vq . 

The collection of spaces V® , j G Z, discussed in Theorem [3] is called 
the quincunx Haar multiresolution analysis in L 2 (Q"), the function (ft from 
axiom (e) is its scaling function. It follows from (J7J) that eft = 1 Z 2. 



It was proved in [6] that the function 

i> = HA-)-t(A--(V / 2 ^ (8) 

is a wavelet function corresponding to the quincunx MRA, i.e. the system 
— a), a € Ip} is an orthonormal basis for the wavelet space Wq := 
V Q . It follows that the system 

{^-a), jeZ,ae/ 2 2 } (9) 

is an orthonormal basis for L 2 ^ 2 .). 

Next we shall show that this basis corresponds also to the separable 
MRA in some sense. Let {Vj}j^z be a separable MRA with p = d = 2, 
tpQ) = tp( 2 ) = l Z2 , i.e. both {^ {1) } jeZ and {vj 2) } j& are Haar MR As. In 
this case, <& = \% % \% 2 = 1 Z 2. Thus, = 0, which yields Vq = V® . Since 

V 1/2 ] 

we have 

—>-((?::))-((?::))• 

It follows that 

{(j){A 2 x-a), a G if} = {0(| - a), a G if }. 

So, W Q = T/ 2 Q = Fx, which yields that W = © , the 
functions ip(- — a),ip(A ■ —a), a G i|, form an orthonormal basis for Wo 
and the dyadic dilations of these functions form an orthonormal basis for 
L 2 (Ql). These relations may be treated as follows: the quincunx Haar MRA 
is the "square root" of the separable Haar MRA. 
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4 Description of wavelet functions 

As was said above, is an orthonormal basis for L 2 (Q|) generated by a 
single wavelet function ip. In contrast to the real setting, p-adic wavelet 
function corresponding to an MRA is not unique. Now we are going to find 
other wavelet functions corresponding to the quincunx Haar MRA. 
First we shall prove some auxiliary facts. 

Lemma 4. Let tp, if) be compactly supported functions with their supports 
contained in B s , s G N, and such that 

$(x) = h a ip(x - a). (10) 
ae/f 

Then h a = whenever | a, 1 2 > 2 s . 

Proof. We rewrite (|10p in the form 

i>{x)= h a ip(x-a) + h a ip(x-a). 

\a\ 2 <2" M 2 >2 S 

The second sum on the right-hand side vanishes whenever \x\2 < 2 s because 
supp^ C B s . The first sum on the right-hand side and the left-hand side 
vanish whenever \x\% > 2 s because supp-0 C B s and suppV' C B s respec- 
tively. So, the second sum is equal to zero for all x 6 Q§. 

Lemma 5. Let n £ Zg, and let ^ be the function defined in ([8]). If An G Z|, 
then + n) = ^, if An Z|, then + n) = — '0. 

Proof. Observe that An G Z 2 if and only if both coordinates of n are 
either even or odd. Therefore, if An £ Z|, then both coordinates of An are 
integers, otherwise both coordinates of An are semi- integers. It remains to 
note that the function <j) is 1 — periodic with respect to each variable. <) 
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Theorem 6. All compactly supported wavelet functions corresponding to the 
quincunx Haar MRA are given by 

2 3 -l 2 S -1 



k=0 1=0 ^ 

where s 6 N, i/j is the function defined by $8$, 

„au 2 S -1 2 S -1 



k/2> 
I /2 s 



OLM 



2S lqp, 



2- 2s {-l) k e~W E E e" 

p=0 <?=0 



p=0 g=0 

7g P G C, |7 9P | = 1, p,q = 0, . . . ,2 s - 1. 



if 1 = 0, 

Iqp, if 



fill 



(12) 



Proof. Let t/j be a compactly supported wavelet function corresponding 
to the quincunx Haar MRA, i.e., the system {ip(- — a), a G /|} is an 
orthonormal basis for Wq , and supp-^ C B s . Since {ip(- — a), a G if} i s a 
basis for W,^ and V> G Wq , we have 

oelf 

Due to Lemma 01 taking into account that supp ip C B\ C B s , supp ip C B s , 
we obtain ([lip with some coefficients ajy. 

Since — a), a G /|} is an orthogonal system, evidently, ^ is orthog- 
onal to — a) whenever a G l\ and a ^ (^j^s ) , A:, / = 0, 1, . . . , 2 s — 1. Thus, 
the system {ip(- — a), a G if} is an orthonormal system if and only if the 
system consisting of the functions ip (^x — (^O) ) = 0, 1, . . . , 2 s — 1, is 
orthonormal. 
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Next we need the following notations. Set 





1 A O . 


. Q Q \ 




( o o . 


. Q -I \ 




O A . 


. o o 




I o . 


. o o 


A = 


O O . 


. A O 


, n = 


O I . 


. o o 




\ o o . 


. OA/ 




\ o o . 


.10/ 



where O, I, A are 2 s x 2 s matrices, O is the zero matrix, I is the unity matrix, 

/ ... -1 \ 



A 



1 
1 



\ 



.. 
.. 



/ 



Let a be a 2 2s -dimensional vector (column) whose iV-th coordinate, N = 
2 s l + k, k, I = 0, 1 ... 2 s — 1, is Ofc/, and let ^, ^ be 2 2s -dimensional vector- 
functions (columns) whose iV-th coordinates, N = 2 s l+k, k, I = 0, 1, . . . , 2 s — 
1, are ip (x — (^0) anc ^ ^ ( x — (t/2- 3 )) res P ec tively. Using Lemma [5l we 
obtain ^ = ', where D is a 2 2s x 2 2s matrix whose iV-th row, N = 2 s l + k, 
k,l = 0,1,... ,2 s - 1, is n l A k a. 

Due to orthonormality of the system {ip( - — a), a £ /|}, the coordinates 
of ^ form an orthonormal system if and only if the matrix D is unitary. 
To describe all unitary matrices D we need to find all a € C 22s such that 
the vectors a, Aa, . . . , £l l A k a, . . . , 2S_1 A 2 " _1 a are orthonormal. We have 
already one such vector ao = (1,0, .. . , 0) T because the matrix Dq whose 
rows are ao, Aao, • • • , Q l A k ao, . . . , Q 2S ~ 1 A 2S ~ 1 ao is, evidently, unitary. 

Let us prove that the rows of D are orthonormal if and only if a = BctQ, 
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where B is a unitary matrix such that 

AB = BA, (13) 

VLB = BQ. (14) 

Indeed, let a = Bao, with B being a unitary matrix satisfying f)13[) . (|14p . In 
this case J7 z A fc a = £l l A k Bao = B(£l l A k ao), and the unitarity of Dq implies 
the unitarity of D. Conversly, if the rows of D are orthonormal, taking into 
account that the rows of -Do are also orthonormal, we conclude that there 
exists a unitary matrix B such that 

D T = BDq (15) 

So we have, in particular, 

a = Ba Q , Aa = BAa , A 2 a = BA 2 a , . . . , A T ~ l a = BA 2 ^ 1 ^. 

Substituting the first equality into the second one, we obtain (AB — BA)ao = 
0. Using this and substituting the first equality into the third one, we obtain 
(AB — BA)(Aao) = 0. After similar manipulations with other equalities we 
have 

(AB - BA)(A p a ) = 0, p = 0, 1, . . . , 2 s - 2. (16) 

Since A 2 " a = —a, A 2 " ao = — a>o, substituting a = Bao and using (fTo]) we 
obtain (AB - BA)(A 2S - 1 a ) = 0. Next it follows from ([15]) that 

na = Ba , UAa = BnAa , VtA 2 a = M7A 2 a , . . . , SlA^^a = BVtA^^ao- 

Substituting a = Bao hito the first equality, we have (QB — BQ)ao = 0. 
Using this, taking into account that QA = Af2 and substituting a = Bao 
into the second equality, we obtain that (AB — BA)({lao) = 0. After similar 
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manipulations with the other equalities and with the equalities £1A 2S a = 
—Qa, QA 2 " ao = —ttao, we have 



(AB - BA)(A p n q a ) = 0, p = 0, 1, . . . , 2 s - 1. 



Continuing this process, we get that (AB — BA)(A p ao) = for all p,q = 
0, 1, . . . , 2 s - 1. Since {A p n%o}p* g =o is a basis for C 22s , we obtain $T5§. In 
a similar way we can check that (|14|) holds. 

Thus we should describe all unitary matrices B satisfying (fl~3|h (fT4"|) . Let 



B 



B, 



oo 



B t 



0,2 S -1 



\ -B2 S -1, ••• ^2 S -1,2 S -1 / 



where Bij is a 2 s x 2 s matrix. It follows from ()14[) that 

' — -S2 S -1,0 — £?2 S -1,1 

-Boo B01 

\ -82^-2,0 82 s -2,1 
This yields 



— B2 S -1,2 S 
Bo,2 s -l 

i?2 s -2,2 s - 



-1 \ / B01 

fill 



Bo,2 s -i 

Bl O'-l 



-Boo \ 
-Bio 



Boo — Bn — 
Bqi = B il+ i 



B2 S -1,2 S -1, 

= B 2 s-l-l t 2 s -l 

-B2"-lfl = ■ ■ ■ = 



'2 S -1J-1 



, V . . . , . 



i.e. 



B 



1 Po Pi 

— /?2 S -1 Po 

V -Pi -Pi 



/?2 s -2 

A) / 
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where = Bqi . It is not difficult to see that any such matrix B satisfies ([T 
Since 



AB 



I A/3 
A(-^-i) 



32 
32 s -2 



BA 



( A) A /3iA ... /3 2 »-iA \ 

(-/3 2 *-i)A /3 A ... ft s _ 2 A 



\ A(-/3i) A(-/3 2 ) ... Ajflb / V (-ySi)A (-ft)A ... /3 A, / 

(PD holds if and only if /3„A = A/3„, for all i/ = 0, . . . , T - 1. It follows 
from the proof of Theorem 4.1 in [17J that /3 n commutes with A if and only 
if f3 v = Cf3 v C , where f3 v is a diagonal matrix, C is a matrix whose entries 
are given by 



c pq = 2- s '\-iy e --^ , p, g = 0, . . . , 2 s - 1. 



Since C is a unitary matrix, £> is unitary if and only if the matrix 



( fa fa ■■■ fc-i \ 



(17) 



B 



2 S -1 



A) • ••/?: 



2 s -2 



V -h 



fa ) 



is unitary. Set Q\ = (A} u; , . . . , A,^ ± ') 1 , where Xt"* is the l-th diagonal 
element of f3 v . It is easy to see that B is unitary if and only if the vectors 
9i,A6i, . . . , A 2S ~ 1 9i form an orthonormal system for all / = 0, . . . , 2 s — 1. It 
follows from the proof of Theorem 4. 1 in [T7| that all such Q\ are given by 

I Ho ... \ 
7/1 ••• 



C 



\ 



C- l e lt 



7/,2 s -i / 
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where j [k € C, \^ k \ = 1, e 1 = (1,0, ...,0) T € M 2 ". Substituting (JT7|), we 
have 



A 



(") 



2 S -1 

2 _ s( _ 1)1/e _^ j2 e-*»^7lfc, /,^ = 0,...,2 S -1. (18) 

fc=0 

Again using (|17p . we obtain 

/c o ... o\ / c- 1 o ... o \ 



C 



B 



C J \ 
( C 



c- 



c 



a 

W io \ 

<^2 S -1 



. c j v y 



, (19) 



where 



/ is N 



V 1 / 



2 -s/2 



\ \ {V) I 
\ A 2 S -1 / 



(20) 



To prove (|12[) it remains to combine (|19p with (|20p . (|18p . and take into 
account that «ofc is the fc-th coordinate of C5q and is the fc-th coordinate 
of C5 2 s-i for Z = 1,...,2 S -1. 
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